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FOREWORD 


This  status  report  of  the  Optical  Sciences  Center,  University  of  Arizona, 
Tucson,  Arizona,  under  U.S.  Air  Force  Contract  FO-469S-67-C-0197,  covers 
research  progress  for  the  poriod  1 January  1975  to  31  March  1975. 

This  contract  was  issued  by  the  Department  of  the  Air  Force,  Space  and 
Missile  Systems  Organization  (AFSC),  P.0.  Box  92960,  Worldway  Postal 
Center,  Los  Angeles,  California  90009,  and  is  administered  by  the  ONRRR, 
University  of  Arizona,  Tucson,  Arizona  85721. 

The  use  of  trade  names  in  this  report  does  not  constitute  an  official 
endorsement  or  approval  of  the  use  of  such  commercial  hardware  or  soft- 
ware. This  report  may  not  be  cited  for  purposes  of  advertisement,  and 
its  publication  does  not  constitute  Air  Force  approval  of  findings  or 
conclusions. 
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PROFESSIONAL  PUBLICATIONS  AND  ACTIVITIES 
OF  OPTICAL  SCIENCES  CENTER  PERSONNEL 
(listed  chronologically) 


ft 
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Journal  Publications 


George  Lamb,  Jr.,  "BSrklund  transformation  for  certain  nonlinear  evolution 
equations,"  J.  Math  Phys.  15(12) : 2157-2165,  Dec.  1974  (publication 
not  received  until  this  quarter). 

D.  Rogovin  and  Marian  0.  Saully,  "Does  the  'two-level  atom'  picture  of 
a Josephson  junction  have  a theoretical  foundation  in  BCS?"  Ann. 
Phys.  88(2) :371-396,  Dec.  1974  (publication  not  received  until  this 
quarter) . 

R.  L.  Shoemaker,  S.  Stenholm,  and  Richard  G.  Brewer,  "Collision- induced 

optical  double  resonance,"  Phys.  Rev.  10(6) : 2037-2050,  Dec.  1974 
(publication  not  received  until  this  quarter). 

S.  F.  Jacobs,  M.  Sargent  III,  J.  F.  Scott,  and  M.  0.  Saully,  eds..  The 

Phyaiae  of  Quantum  Electronics , Vol.  II,  Laser  Applications  to  Optics 
arid  Spectroscopy,  Addison-Wesley,  Reading,  Mass.,  1974  (publication 
not  received  until  this  quarter). 

Peter  H.  Bartels,  George  B.  Olson,  Jack  M.  Layton,  Robert  E.  Anderson, 
and  George  L.  Wied,  "Computer  discrimination  of  T and  B lympho- 
cytes," Acta  Cytologica  19(1) : 53-57 , Jan. -Feb.  1975. 

R.  E.  Anderson,  G.  B.  Olson,  C.  Shank,  T.  L.  Howarth,  G.  L.  Wied,  and 
P.  H.  Bartels,  "Computer  analysis  of  defined  populations  of  lympho- 
cytes irradiated  in  vitro.  I.  Evaluation  of  murine  thoracic 
duct  lymphocytes,"  Acta  Cytologica  19(1) :126-135,  Jan. -Feb.  1975. 

W.  M.  Fairbarik,  Jr.,  G.  K.  Klauminzer,  and  A.  L.  Schawlow,  "Excited-state 
absorption  in  ruby,  emerald,  and  MgO:Cr3+,"  Phys.  Rev.  Bll:60-76, 

Jan.  1975. 

V.  .V.  Mahajan  and  J.  D.  Gaskill,  "Diffraction  of  light  by  sound  waves  in 
optically  isotropic  media,"  Opt.  Eng.  14:91-97,  Jan. -Feb.  1975. 

A.  S.  Marathay,  "Polarized  light  and  its  application,  I,"  Opt.  Eng. 

14 : S- 17  to  S-21 , Jan. -Feb.  1975. 

M.  P.  Rimmer  and  J.  C.  Wyant , "Evaluation  of  large  aberrations  using  a 
lateral-shear  interferometer  having  variable  shear,"  Appl.  Opt. 
14:142-150,  Jan.  1975. 
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W.  M.  Fairbank , Jr.,  T.  W.  ilUnsch,  and  A.  L.  Schawlow,  "Absolute  mea- 
surement of  very  low  sodium-vapor  densities  using  laser  resonance 
fluorescence,"  J.  Opt.  Soc.  Am.  65(2) : 199-204,  Feb.  1975. 

J.  B.  Hambenne  and  M.  Sargent  III,  "Physical  interpretation  of  bistable 
unidirectional  ring-laser  operation,"  IEEE  J.  Quantum  Electronics 
QE-ll:90-92,  Feb.  1975. 

D.  Rogovin,  H.  Tigalaar,  and  M.  Sargent  III , "A  simple,  reliable  moael 
of  the  bending  mode  of  water,"  Chem.  Phys.  Letters  31:147-153,  Feb. 
1975. 

H.  Shih,  M.  0.  Scully,  P.  V.  Avizonis,  and  W.  H.  Louisell,  "Multimode 
approach  to  the  physics  of  unstable  laser  resonators,"  Phys.  Rev. 
All:630-647,  Feb.  1975. 

0.  N.  Stavroudia , Book  review  of  Optics  by  E.  Hecht  and  A.  Zajac, 

Addison-Wesley,  Reading,  Mass.,  1974;  Phys,  Today  28(2): 54-55,  Feb. 
1975. 

F.  0.  Bartell  and  W.  L.  Wolfe,  "New  approach  for  the  design  of  blackbody 
simulators,"  Appl.  Opt.  14(2) : 249-252,  Feb.  1975. 

R.  A.  Johnston  and  C.  L.  Babcock,  "Composition  dependence  of  elastic 
moduli  in  Na20-Ti02-Si02  glasses,"  J.  Am.  Ceram.  Soc.  58(3-4): 

85-87,  Mar. -Apr.  1975. 

V.  N.  Mahajan,  "Real-time  wavefront  correction  through  Bragg  diffraction 
of  light  by  sound  waves,"  J.  Opt.  Soc.  Am.  65(3) : 271-278,  Mar.  1975. 

D.  Rogovin,  "Electrodynamics  of  Josephson  junctions,"  Phys.  Rev.  Bll(5): 
1906-1908,  Mar.  1975. 

M.  H.  Farmelant,  C.  C,  DeMeester,  D.  T.  Wilson,  and  H.  H.  Barrett,  "Ini- 
tial clinical  experiences  with  a Fresnel  zone  plate  imager,"  J. 

Nucl.  Med.  16:183-187,  1975. 

D.  Rogovin,  "Phase  coherence  in  two-band  superconductors,"  J.  Low 
Temp.  Phys.  18(i):l-2,  1975. 
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Invited  Papers 

B.  R.  Frleden,  "Super-resolution,  entropy  and  phase  contrast,"  at  a 

special  seminar  of  the  National  Radio  Astron.  Observ.,  Charlottes- 
ville, Virginia,  Jan.  13,  1975. 

B.  0.  Seraphin,  "Solar  utilization  now  short  course  (SUN)  on  solar 

energy,"  Arizona  State  University,  Tempe,  Arizona,  Jan.  13-14,  1975. 

M.  0.  Scully,  "Charge  exchange  excitation--x-ray  laser,"  Orbis 

Scientiae  II,  Center  for  Theoretical  Studies,  University  of  Miami, 
Coral  Gables,  Florida,  Jan.  19-24,  1975. 

M.  0.  Scully,  "Resonant  charge  exchange,"  Orbis  Scientiae  II,  Center 

for  Theoretical  Studies,  University  of  Miami,  Coral  Gables,  Florida, 
Jan.  19-24,  1975. 

B.  0.  Seraphin,  "Spectrally  selective  surfaces  for  high  temperature 
photo  thermal  conversion,"  1975  Joint  Meeting  of  the  American 
Physical  Society  and  the  American  Association  of  Physics  Teachers, 
Anaheim,  California,  Jan.  29-Feb.  1,  1975. 

M.  0.  Scully,  "Resonant  charge  exchange,"  and  "Short  X and  superradiant 
resonators,"  5th  Winter  Colloquium  on  Quantum  Electronics,  Snowmass, 
Colorado,  Feb.  1-6,  1.175. 

W.  Swindell,  L,  R.  Doose,  M.  Tomasko,  and  J.  Fountain,  "The  Pioneer  11 
images  of  Jupiter,"  presented  at  the  6th  Annual  Meeting  of  the 
Division  for  Planetary  Sciences,  American  Astronomical  Society, 

Feb.  17-21,  1975. 
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Conference  Proceedings 


H.  H.  Barrett,  "Fresnel  zone  plate  and  related  coded  apertures,"  in 
Proceedings  of  the  SPIE  National  Seminar  cn  Coherent  Optical 
Processing,  Vol.  52,  H.  J.  Caulfield,  ed.  (Meeting  held  August 
21-22,  1974,  San  Diego,  California). 

R.  Buchroeder,  "Tilted  component  optical  design  methods  for  infrared."  in 
Proceedings  of  the  SPIE  National  Seminar  on  Applications  of  Geo- 
metrical Optica,  Vol.  39,  Warren  J.  Smith,  ed.  (Meeting  held  August 
27-29,  1973,  San  Diego,  California). 

James  J.  Burke,  "Estimating  objects  from  their  blurred. and  grainy 

images,"  in  Proceedings  of  the  Technical  Program,  Electro-Optical 
Systems  Design  Conference-1974  West,  International  Laser  Exposition 
(Meeting  held  November  5-7,  1974,  San  Francisco,  California). 

A.  B.  Meinel,  "Optical  interfaces  in  solar  energy  utilization,"  in 

Proceedings  of  the  SPIE  National  Seminar  on  Effective  Systems  In- 
tegration & Optical  Design,  Vol.  54,  G.  Wilkerson  and  R.  W.  Poindexter, 
eds.  (Meeting  held  August  21-23,  1974,  San  Diego,  California). 

M.  Ruda,  "Optical  tolerancing  for  multi-mirror  telescope,"  in  Pro- 
ceedings of  the  SPIE  National  Seminar  on  Applications  of  Geometrical 
Optics,  Vol.  39,  Warren  J.  Smith,  ed.  (Meeting  held  August  27-29, 

1973,  San  Diego,  California). 

R.  A.  Schowengerdt,  "Measurement  of  earth  resources  technology  satellite 
(ERTS-1)  multispectral  scanner  OTF  from  operational  imagery,"  in 
Proceedings  of  the  SPIE/OSA  National  Seminar  on  Image  Assessment  & 
Specification,  Vol.  46,  D,  Dutton,  ed.  (Meeting  held  May  20-22,  1974, 
Rochester,  New  York) . 

R.  V.  Shack,  "Significance  of  phase  transfer  function,"  in  Proceedings 

of  the  SPIE/OSA  National  Seminar  on  Image  Assessment  & Specification, 
Vol.  46,  D.  Dutton,  ed.  (Meeting  held  May  20-22,  1974,  Rochester, 

New  York) . 

R.  V.  Shack,  "Use  of  normalization  in  application  of  simple  optical 

systems,"  in  Proceedings  of  the  SPIE  National  Seminar  on  Effective 
Systems  Integration  & Optical  Design,  Vol.  54,  G.  Wilkerson  and  R.  W. 
Poindexter,  eds.  (Meeting  held  August  21-23,  1974,  San  Diego, 
California) . 


R.  V.  Shack,  "Analytic  system  design  with  pencil  5 ruler--advantages  of 
y-y  diagram,"  in  Proceedings  of  the  SPIE  National  Seminar  on  Ap- 
plications of  Geometrical  Optics,  Vol.  39,  Warren  J.  Smith,  ed. 
(Meeting  held  August  27-29,  1973,  San  Diego,  California). 

R.  R.  Shannon,  "Fringe  benefits  in  testing,"  in  Proceedings  of  the  SPIE 
National  Seminar  on  Applications  of  Geometrical  Optios,  Vol.  39, 

Warren  J.  Smith,  ed.  (Meeting  held  August  27-29,  1973,  San  Diego, 
California) . 

k.  R.  Shannon,  "Massive  optical  technology  at  the  Optical  Sciences 

Center,"  in  Proceedings  of  the  Technical  Program,  Electro-Optical 
Systems  Design  Conference- 19 74  West,  International  Laser  Exposition 
(Meeting  held  November  5-7,  1974,  San  Francisco,  California). 

Philip  N.  Slater  and  William  L.  Wolfe,  eds..  Proceedings  of  the  SPIE 
National  Seminar  on  Scanners  & Imagery  Systems  for  Earth  Observa- 
tion, Vol.  51,  (Meeting  held  August  19-20,  1974,  San  Diego, 

California) . 

P.  N.  Slater,  "Specifications  for  remote  sensing  photographic  systems," 
in  Proceedings  of  the  SPIE  National  Seminar  on  Effective  Systems 
Integration  & Optical  Design,  Vol.  54,  G.  Wilkerson  and  R.  W. 
Poindexter,  eds.  (Meeting  held  August  21-23,  1974,  San  Diego, 
Crlifomia) . 

James  4.  Palmer  and  William  L.  Wolfe,  "The  optics  and  calibration  of 
tl  e Venus  solar  flux  radiometer,"  in  Proceedings  of  the  Technical 
Pogrom,  Electro-Optical  Systems  Design  Conference-1974  West, 
International  Laser  Exposition  (Meeting  held  November  5-7,  1974, 

San  Francisco,  California). 

Frederick  0.  Bartell  and  William  L.  Wolfe,  "A  comparison  of  the  performance 
of  blackbody  simulators  of  different  shapes,"  in  Proceedings  of 
the  Technical  Program,  Electro-Optical  Systems  Design  Conference- 
1974  West,  International  Laser  Exposition  (Meeting  held  November 
5-7,  1974,  San  Francisco,  California). 
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Contributed  Papers 


R.  E.  Hahn  and  B.  0.  Seraphin,  "Thick  semiconductor  films  for  ohoto- 

thermal  solar  energy  conversion,"  Conference  on  Structure  Property 
Relationships  in  Thick  Film  and  Bulk  Coatings,  San  Francisco, 
California,  Feb.  10-12,  1975. 

J.  Hambenne  and  M.  Sargent  III,  "Physical  interpretation  of  bistable, 
unidirectional  ring  laser  operation,"  Bull.  Am.  Phvs.  Soc., 
Anaheim,  California,  1975. 


Technical  Reports  and  Government  Reports 

0.  Seraphin,  "Chemical  vapor  deposition  research  for  fabrication  of 
solar  energy  convertors,"  NSF/RANN  1974  Annual  nrogress  renort, 

(•rant  No.  GI-36731X. 

A.  Schowengerdt  and  P.  N.  Slater,  "Evaluation  of  ERTS-1  image  sensor 
spatial  resolution  in  photographic  form,"  NASA/Goddard  Space 
Flight  Center,  Jan.  1975. 

A.  Franken,  "Semi-annual  FY7S  applied  optics  research  report," 

SAMSO  TR  75-112,  January  1975. 

R.  Shannon,  R.  V.  Shack,  T.  R.  Gurski,  D.  W,  Hillman,  F.  A.  Honf, 

R.  L.  Shoemaker,  I.  R.  Harrison,  and  J.  C.  Wyant,  "Optical  tracking 
system  study  interim  report,"  Space  and  Missile  Test  Center, 
FO-4703-75-C-0210,  February  1975. 

H.  Bartels,  J.  J.  Burke,  B.  R.  Frieden,  R.  S.  Hershel,  and  L.  Whoeler, 
"Basic  and  applied  research  in  optical  science  and  technology," 
annual  report  on  F-33657-C-605,  February,  1975. 
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Colloquia 

l).  N.  Rogovin,  "Molecular  spectroscopy,"  Kirtland  AFB,  January  8-9, 

1975. 

M.  Sargent  III,  "Multimode  model  of  water  molecule,"  Kirtland  AFB, 
January  8-9,  1975. 

G.  Lamb,  Jr. , "Coherent  optical  pulse  propagation,"  Applied  Mathematics 

Seminar,  California  Institute  of  Technology,  January  20,  1975. 

B.  R.  Fricden,  "Super-resolution  by  Fresnel-encoding  the  image," 

Optical  Sciences  Center  Colloquium,  January  30,  1975. 

H.  II.  Barrett,  "Coded  aperture  imaging,"  Biomedical  Engineering  Seminar, 

University  of  Arizona,  February  3,  1975. 

R.  I..  Shoemaker,  "Infrared  quantum  beat  spectroscopy,"  Physics  Dent., 
University  of  Arizona,  February  18,  1975. 

.).  D.  Gaskill,  "The  status  and  future  of  holography,"  Optical  Sciences 
Center  Colloquium,  February  20,  1975. 

M.  0.  Scully,  invited  lectures  at  Physics  Dept.,  University  of  Colorado 
end  Joint  Institute  for  Laboratory  Astrophysics,  Feb.  24-28,  1975. 

M.  Meinel,  "Solar  energy,"  University  of  Wisconsin,  Milwaukee,  Wisconsin 
March  5,  1975. 

W.  Fairbank,  "Detection  of  small  quantities  of  matter  by  laser  resonance 
fluorescence,"  Colorado  State  University,  March  11,  1975. 

.1.  M.  Palmer  and  M.  G.  Tomasko,  "The  1978  Pioneer-Venus  solar  flux 

experiment,"  University  of  Arizona  Atmospheric  Sciences  Colloquium, 
March  14,  1975. 

A.  B.  Meinel  and  M.  Meinel,  "Solar  energy,"  Weber  State  College,  Ogden, 
Utah,  March  27,  1975. 

I . Wheeler,  "Review  of  psychophysical  investigations  of  image  quality," 

SAMSO,  Los  Angeles,  California,  March  12,  1975. 
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Organization  of  Scientific  Meetings  and  Courses 

P.  A.  Franken,  Chairman,  "Committee  on  Society  Objectives  and  Policy 
(SOAP)  of  the  Optical  Society  of  America,"  Tucson,  Arizona 
Jan.  10-11,  1975. 

M.  0.  Scully,  Session  moderator,  "X-ray  laser  session,  Orbis  Scientiae 
II,"  at  Center  for  Theoretical  Studies,  University  of  Miami, 

Coral  Gables,  Florida,  Jan.  19-24,  1975. 

B.  0.  Seraphin,  organized  session  on  "Physics  of  solar  energy  conversion," 
at  American  Physical  Society  meeting,  Anaheim,  California,  Jan. 

29- Feb.  1,  1975. 

M.  0.  Scully,  Co-director,  "5th  Winter  Colloquium  on  Quantum  Electronics," 
Snowmass,  Colorado,  Feb.  1-6,  1975. 

B.  0.  Seraphin,  organized  planning  session  of  the  "National  Science 

Foundation's  Workshop  on  Material  Sciences,"  at  Harvard  University, 
Cambridge,  Massachusetts,  Feb.  13-17,  1975. 

R.  V.  Shack,  "Optical  Fabrication  fi  Testing  Workshop,"  sponsored  by 
the  Optical  Society  of  America,  San  Francisco,  California,  Mar. 

13-16  and  Anaheim,  California,  Mar.  19-21,  1975. 

R.  R.  Shannon,  Chairman,  Panel  discussion  on  "Applications  of  Laser 
Spectroscopy,"  at  Optical  Society  of  America  meeting,  Anaheim, 
California,  March  20,  1975. 
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New  Committee  Memberships 

P.  A.  Franken,  Ballistic  Missile  Defense  Technology  Advisory  Panel, 
Department  of  the  Army,  Huntsville,  Alabama. 

M.  Sargent  III,  U.S.  Senior  Scientist  Award,  fellowship  for  pursuit  of 
research  in  Germany. 

B.  0.  Seraphin,  participant  in  joint  USSR  - USA  workshop  on  Solar 
Central  Power  Stations,  Colorado  Stale  University,  Ft.  Collins, 
Colorado,  March  30-April  1,  1975. 

B.  0.  Seraphin,  Consultant  to  NASA  - Marshall  Space  Flight  Center  on 
recommendations  for  the  implementation  of  the  Solar  Heating  § 
Cooling  Act  being  put  together,  March  1975. 

L.  Wheeler,  Acting  Chairman,  Department  of  Psychology,  University  of 
Arizona,  January  1,  1975. 

W.  L.  Wolfe,  American  Editor  of  Infrared  Physics. 

W.  L.  Wolfe,  Chairman  of  Optical  Society  of  America  Technical  Group  on 
Radiometry . 

J.  C.  Wyant,  Guest  editor  of  Sept. -Oct.  issue  of  Optical  Engineering, 
organi'ing  the  special  issue  on  "Application  of  Holography." 
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Public  Service  Speaking 

W.  Swindell,  "Pioneer  10  and  ll--first  pictures  from  Jupiter,"  Optical 
Sciences  Center  Public  Evening,  Jan.  13,  1975. 

B.  0.  Seiaphin,  "Spectrally  selective  coatings  for  solar  energy  con- 
version," Tucson  Section  Optical  Society  of  America,  January  27,  1975. 

B.  0.  Seraphin,  participat  d in  a press  conference  on  "Solar  energy 
conversion,"  Anaheim,  California,  January  30,  1975. 
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RESEARCH  ACTIVITIES 


IMAGE  PROCESSING 

Digital  Processing  of  Annular  Coded-Aperture  Imagery 
(R.  G.  Simpson,  H.  H.  Barrett,  and  H.  D.  Fisher) 


Summary 


In  the  past,  annular  coded-aperture  images  have  been  reconstructed 
by  correlating  the  annulus  with  the  coded  image.  The  basic  improvement 
suggested  in  this  report  is  the  addition  of  a linear  radial  frequency 
weighting  in  the  Fourier  plane.  Reconstruction  of  point  and  disk 
objects  was  simulated  with  a computer  program.  The  results  show  the  ad- 
vantage of  this  modification  in  the  processing  scheme.  When  the  assumption 
is  made  that  the  detector  is  an  Anger  camera,  the  resolution  obtained 
with  the  improved  processing  of  the  coded  image  is  equal  to  that  obtained 
with  conventional  apertures. 

The  annulus  is  intermediate  between  the  pinhole  and  the  Fresnel 
zone  plate  with  regard  to  both  collection  efficiency  and  the  number  of 
counts  required  for  a given  signal-to-noise  ratio  (SNR).  It  therefore 
offers  an  improvement  over  pinhole  apertures  without  demanding  the  in- 
creased count  rate  and  resolution  required  of  detectors  when  Fresnel 
zone  plate  coded  apertures  are  used. 

Introduction 

The  imaging  problem  encountered  in  nuclear  medicine  is  the 
formation  of  a picture  of  the  distribution  of  radiopharmaceuticals  within 
a patient.  Isotopes  commonly  used  in  nuclear  medicine  emit  gamma  rays 
with  energy  in  the  100  to  400  keV  range.  For  photons  in  this  energy 
range,  reflection,  refraction,  and  diffraction  effects  are  so  small  that 
they  are  practically  useless  for  imaging  purposes.  The  only  alternative 
is  to  use  the  absorption  of  gamma  rays  by  matter  to  cast  shadows. 
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One  simple  shadow-casting  device  is  the  pinhole.  Another  is  the 
collimator,  which  restricts  the  field  of  view  of  a given  section  of  the 
detector  to  a given  region  of  the  object.  When  either  of  these  apertures 
is  used,  a compromise  must  be  made  between  collection  efficiency  and 
resolution;  as  one  improves,  the  other  necessarily  degrades. 

Coded  apertures  of  various  descriptions  have  been  proposed  as 
alternatives  to  the  pinhole  and  the  collimator.  While  coded  apertures 
collect  far  more  counts  than  conventional  apertures,  they  also  require 
more  counts  to  obtain  the  same  SNR.  One  measure  of  a potential  coded 
aperture  is  whether  the  increase  in  collected  counts  exceeds  the  required 
increase  in  counts  for  a constant  SNR.  In  general,  there  will  be  a net 
advantage  to  coded  apertures  for  relatively  small  objects,  but  the 
advantage  will  disappear  as  the  object  size  approaches  the  full  field 
of  view. 

One  class  of  coded  apertures,  referred  to  as  filled  apertures, 
is  roughly  50%  transparent  to  the  gamma  rays.  The  Fresnel  zone  plate  and 
Dicke's  random  pinhole  array  are  examples  of  this  type  of  aperture. 

Another  class  of  coded  apertures,  called  dilute  apertures,  is  much  less 
than  50%  transparent.  The  nonredundant  pinhole  array  is  one  example  of 
this  class.  Another  is  the  annulus  that  will  be  described  here 
Dilute  apertures  are  intermediate  between  the  filled  apertures  and  the 
pinhole  with  regard  to  both  collection  efficiency  and  the  number  of  counts 
needed  for  a given  SNR. 

One  additional  feature  of  coded  apertures  is  their  tomographic 
capability.  This  refers  to  the  fact  that  different  planes  in  the  image 
can  be  brought  into  focus  separately.  The  scale  of  a shadow  created 
by  an  object  point  will  depend  on  the  distance  of  the  object  from  the 
plane  of  the  aperture.  In  this  way  depth  information  is  stored  in  the 
coded  image. 

While  this  tomographic  capability  is  a useful  property,  it 
causes  some  difficulty.  Reconstruction  usually  involves  correlation  of 
the  coded  image  with  a processing  function.  The  scale  of  the  processing 
function  is  chosen  so  that  a certain  plane  in  the  image  will  be  cor- 


rectly  reconstructed.  However,  shadows  with  many  different  scales  are 
present  in  the  coded  image.  Contributions  to  the  reconstruction  due  to 
these  other  shadows  must  not  add  structure  to  the  image  that  might 
result  in  a faulty  diagnosis.  Effects  due  to  shadows  with  the  incorrect 
scale  will  be  referred  to  as  out-of-focus  behavior. 

The  Annulus 

The  annulus  consists  of  a sheet  of  absorbing  material  such  as 
lead  with  a circular  ring  cut  out  of  it.  The  shadow  created  by  such 
an  aperture  will  be  described  by  an  inner  radius  rj  and  an  outer 
radius  r2. 

Walton  first  used  the  annulus  as  a coded  aperture  in  nuclear 
medicine.1  Other  investigators  have  taken  an  interest  in  it  since 
then.2*3  Walton's  reconstruction  process  was  effectively  a correlation 
of  the  coded  image  with  a t lin  annulus.  The  point  spread  function  (PSF) 
associated  with  such  an  encoding-reconstruction  scheme  is  shown  in  Fig.  1. 
The  undesirable  feature  of  this  function  is  its  rather  long  tail.  Object 
points  are  then  capable  of  influencing  the  reconstruction  at  points 
rather  far  from  their  actual  positions.  For  this  reason  it  is  desirable 
to  have  the  PSF  as  compact  as  possible.  One  way  to  make  the  annular 
aperture  PSF  more  compact  will  now  be  discussed. 


>• 


Fig.  1.  Point  spread  function  corresponding 

to  the  autocorrelation  of  the  annulus. 
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If  one  approximates  the  annulus  by  a ring  delta  function, 

6(r  - /?) , the  autocorrelation  has  the  form4 

AT) 2 

6(r  - R)  * * 6(r  - R)  = — — — - . (1) 

r(4 R2  - r2)'5 

(The  double  star  notation  indicates  a two-dimensional  correlation 
operation.)  Away  from  the  singularity  at  r = 2R , the  function  behaves 
as  I//*.  It  is  this  behavior  that  needs  to  be  eliminated.  A trip  to 
the  frequency  domain  provides  the  key  to  this  problem. 

The  Fourier  transform  of  a thin  annulus  is  closely  approximated 
by  2nr*</0  (2itr*p)  where  ? = (rj  * r2)/2  and  p is  the  magnitude  of  the  spatial 
frequency  vector  in  cycles  (or  line  pairs)  per  unit  length.  To  a 
good  approximation,  beyond  the  first  zero, 

U 

2vvJq  (2irrp)  = 2(r*/p)  2 cos(2irr*p  - m/4) . (2) 

The  Fourier  transform  of  the  autocorrelation  of  a thin  annulus  is  then 

4n2r,2e/Q2(2tn,p)  - — cos2  (2irr*p  - h/4) 

2r  — 

= — [1  + cos(4mnp  - u/2) ] 

= ^ [1  + sin  4Ttrp].  (3) 

If  all  the  frequency  components  are  weighted  by  p,  the  frequency  spectrum 
becomes  2r[l  + sin  4tt?p].  The  constant  term  in  frequency  space  will 
result  in  a spatial  function  sharply  peaked  about  the  origin.  The  sine 
term  will  contribute  to  the  spatial  function  for  values  of  r near  2?. 

That  the  1/r  behavior  of  the  PSF  can  bo  eliminated  in  this  way  is  not 
surprising  when  one  recalls  that5 


F{l/r} 


1/P. 


(4) 
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This  use  of  a linear  frequency  weighting  also  occurs  in  the  literature 
on  computerized  axial  tomography  where  it  is  called  a rho-filter. 

It  is  important  to  realize  that  this  linear  weighting  of  fre- 
t quencies  cannot  go  on  indefinitely.  Noise  considerations  require  a 

termination  of  this  boosting.  The  modulation  transfer  function  (MTF) 
of  the  detector  could  be  used  to  damp  the  linear  weighting  at  high 
spatial  frequencies.  Alternatively,  an  artificial  apodization  function 
t F(p)  could  be  included  in  the  frequency  domain  processing. 

Computer  Simulations 

. The  effect  of  the  rho-filter  has  been  investigated  with  the  use 

* 

of  a computer  program  that  performs  Fourier  transforms  of  circularly 
symmetric  functions.  The  coded  image  g(r ) is  given  by  the  convolution 
of  the  object  o(r)  with  the  annulus  ann(r) 

f 

g{ r)  * o(r)  * * ann(r).  (5) 

(Scale  factors,  as  discussed  in  Ref.  6,  have  been  suppressed  in  this 
equation.)  The  coded  image  is  Fourier  transformed  yielding,  for  cir- 
cularly symmetric  functions, 

C(p)  = 0( p)  • ANN (p ) . (6) 

► This  function  is  then  multiplied  by  the  transform  of  a thin  annulus, 

2tfr«7o  (2irrp) , the  rho-filter,  and  the  apodization  function  F(p) 

C'(p)  * 0(p)  • ANN(p)  • 2t\pJq  (2wrp)  • p • F(p).  (7) 

I 

Then  G' (p)  is  inverse  transformed  to  give  the  reconstructed  object  g'(r). 

For  all"  studies  reported  here,  the  dimensions  of  the  annular 
shadow  in  the  detector  plane  were  r1  = 47.5  mm  and  r2  = 52.5  mm.  The 
PSF  with  an  apodization  F(p)  described  by 
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_ (1  for  0 < p < 0.125  lp/mm  ^ 

(0  for  p > 0.125  lp/mm 

is  shown  in  Fig.  2a.  The  1/r  tail  of  Fig.l  has  been  eliminated.  The 
ringing  is  due  to  the  sharp  cutoff  in  the  frequency  plane. 

A gaussian  apodization  that  approximates  the  MTF  of  an  Anger 
camera  with  8-mm  resolution  was  used  in  Fig.  2b.  The  ringing  has  been 
virtually  eliminated.  However,  the  bipolar  wiggle  at  r = 2?  due  to  the 
sin4itrp  term  in  Eq.  (3)  is  still  present.  Except  for  the  wiggle,  the 
PSF  is  now  much  more  compact  and  even  has  a resolution  as  good  as  that 
of  the  camera. 


Fig.  2.  Point  spread  function  obtained  with 
the  rho-filter  included,  (a)  Sharp 
cutoff  in  frequency  plane,  (b)  Gaussian 
apodization  in  frequency  plane. 


The  improvement  afforded  by  the  rho-filter  operation  is  shown 
remarkably  when  disk  objects  are  imaged.  Figure  3 shows  the  reconstruction 
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obtained  if  the  rho-filter  operation  is  omitted.  The  vertical  scale  is 
in  arbitrary  units  bu*  is  the  same  for  each  of  the  disks.  The  magnitude 
at  the  center  depends  linearly  on  the  radius  of  the  disk.  The  effect  of 
the  rho-filter  is  shown  in  Fig.  4.  The  magnitudes  at  the  centers  of  the 
disks  are  now  all  the  same.  The  true  edge  of  each  disk  occurs  re- 
markably close  to  the  half-peak  value  of  the  reconstruction.  The  data 
beyond  r ■ 50  mm  are  due  to  the  bipolar  wiggle  of  the  PSF. 


Fig.  3.  Reconstruction  of  disk  objects  of  several 

sizes  by  correlation  with  an  annulus  without 
the  rho-filter  operation. 


Fig.  4.  Reconstruction  of  disk  objects  of  several 
sizes  by  correlation  with  an  annulus  with 
the  rho-filter  operation. 
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Again,  assuming  a gaussian  model  for  the  Anger  camera  MTF,  the 
reconstruction  of  three  sets  of  concentric  annuli  was  simulated  with  the 
results  shown  in  Fig.  5.  The  radius  of  the  central  disk  was  chosen  to 
be  three  quarters  the  width  of  the  rings  so  that  a well  defined  radial 
frequency  would  result.  Three  quarters  was  chosen  because  the  first 
zero  of  the  Jq  function  occurs  at  approximately  three  quarters  of  the 
distance  between  consecutive  zeros  of  the  Jq  function.  In  Fig.  Sc, 
the  rings  are  slightly  narrower  than  the  resolution.  The  reconstruction 
shown  in  Fig.  5d  still  has  a contrast  of  35%.  For  the  5-mm-wide  rvngs 
shown  in  Fig.  5e,  the  contrast  in  the  reconstruction  is  only  10%. 


h n n 

a 
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h n n 
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\/\/\  d 

e 

EL 

\aA 

Radius 

Fig.  5-  Concentric  ring  objects  and  their  reconstruc- 
tion. The  rho-filter  is  included  along  with 
the  gaussian  MTF  of  the  Anger  camera.  The 
field  of  view  has  been  limited  so  that  arti- 
facts due  to  the  wiggle  at  2 r are  outside 
the  field,  (a)  Object  with  5-mm  ring  width, 
(b)  Reconstruction  of  (a),  (c)  Object  with 

7-mm  ring  width,  (d)  Reconstruction  of  (c) . 
(ej  Object  with  10-mm  ring  width,  (f)  Recon- 
struction of  (ej . 


As  mentioned  earlier,  for  the  tomographic  capabilities  of  the 
coded- aperture  reconstruction  technique  to  be  useful,  the  out-of- focus 
response  should  not  cause  confusing  structure  to  appear.  A parameter 
called  FOCUS,  defined  as  the  ratio  of  the  reconstructing  annular  radius 
to  the  coding  annular  radius,  was  investigated  for  its  effects  on  the  PSF. 
The  results  are  shown  in  Fig.  6.  The  effect  of  FOCUS  on  the  reconstruction 
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of  a disk  object  is  illustrated 
in  Fig.  7.  At  the  present 
time  there  is  no  hard  and  fast 
criterion  for  judging  the  out-of- 
focus  behavior  of  a coded  aperture. 
These  results  do  not  appear  dis- 
turbing, but  at  the  same  time  they 
are  not  ideal. 

Signal-to-Noise  Ratio 

When  a coded  aperture  is 
used  to  image  a point  source,  more 
counts  are  collected  than  for  a 
pinhole.  The  SNR  is  therefore 
better  for  the  coded  aperture  case. 
However,  when  more  than  one  point 
is  being  imaged,  the  shadows  created 
by  the  object  points  may  overlap. 
Then  counts  due  to  one  object  point 
will  be  a noise  source  for  another 


r c nc_  object  point.  The  effects  of  this 
Fig.  6.  Effects  of  FOCUS  on  PSF. 

noise  will  decrease  the  advantage 
enjoyed  by  coded  apertures.7 

The  magnitude  of  this  decrease  depends  on  the  amount  of  overlap. 
Consider  a Fresnel  zone  plate  with  the  same  resolution  as  the  annulus 
described  above.  When  the  object  consists  of  two  points  whose  shadows 
are  displaced  with  respect  to  each  other  by  half  the  maximum  radius  of 
the  zone  plate,  the  overlap  is  roughly  40%  of  the  open  area.  For  the 
annulus  under  the  same  conditions,  the  overlap  is  only  10%  of  the  open 
area.  The  SNR  for  the  zone  plate  is  therefore  degraded  more  than  for  the 
annulus,  but  the  zone  plate  has  the  higher  geometric  collection  efficiency. 
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Fig.  7.  Effects  of  FOCUS  on  the  reconstruction  of  a 
15-nun-radius  disk. 

The  two  factors  work  in  opposition  to  one  another,  suggesting  that  the 
annulus  is  a good  alternative  to  the  zone  plate  for  intermediate  values 
of  the  ratio  of  object  area  to  resolution  cell  size. 

An  exact  comparison  of  the  SNR  characteristics  of  the  annulus 
with  other  coded  apertures  requires  the  calculation  of  the  noise  kernel. 
A computer  program  is  being  written  to  provide  this  information.  The 
value  of  the  annulus  as  a coded  aperture  will  ultimately  depend  on  this 
comparison. 
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Detector  Requirements 

A short  note  needs  to  be  made  here  regarding  an  advantage  of  the 
annulus  over  the  Fresnel  zone  plate  when  an  Anger  camera  is  used  as  the 
detector.  As  discussed  above,  the  annulus  collects  fewer  counts  than 
th*  Fresnel  zone  plate  but  also  needs  fewer  for  the  same  SNR.  Therefore, 
the  count  rate  capability  of  the  camera  is  not  taxed.  In  addition,  the 
required  resolution  is  not  as  great  as  when  an  off-axis  Fresnel  zone 
plate  is  used.  Therefore,  camera  limitations  are  not  as  severe  for  the 
annulus  as  for  the  zone  plate. 

Conclusions 

The  improvement  in  the  PSF  provided  by  the  rho-filter  operation, 
along  with  the  striking  results  in  the  disk  object  reconstruction, 
make  the  annulus  an  attractive  aperture.  The  annulus  appears  well 
suited  for  the  Anger  camera,  and  resolution  as  good  as  that  of  the 
camera  is  attainable.  Finally,  on  empirical  grounds,  the  annulus  shows 
promise  of  being  as  good  as  the  zone  plate  for  values  of  object  area/ 
resolution  cell  size  commonly  encountered  in  nuclear  medicine. 
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MATHEMATICAL  OPTICS 

Design  Modules  and  Optical  Calculations 
(0.  N.  Stavroudis) 

These  two  research  areas  consist  of  developing  new  concepts  and 
their  realization  in  the  form  of  computer  programs  in  the  area  of 
modular  optical  design.  In  addition,  general-purpose  time-sharing 
computer  programs  are  being  developed  for  optical  analysis  and  optical 
design. 

Optical  Design  Modules 

In  the  past  quarterly  report  we  described  how  modules  became 
concentric  at  a particular  value  of  the  shape  parameter.  This  value  of 
the  shape  parameter  is  now  included  in  our  list  of  critical  values 
calculated  in  our  module  programs.  Concentric  modules  can  be  used  to 
form  concentric  lenses.  Such  lenses  have  the  useful  properties  of 
any  concentric  lens  --  zero  coma  and  astigmatism;  zero  distortion  and 
lateral  color.  In  addition,  third-order  spherical  aberration  is  zero. 
However,  Petzval  and  therefore  field  curvature  tend  to  be  large.  For 
the  camera  lens  configuration  axial  color  can  be  made  equal  to  zero  in 
the  manner  described  in  earlier  reports.  However,  Petzval  cannot.  The 
equation  for  zero  Petzval,  in  this  instance,  turns  out  to  be  insoluable. 

The  reason  for  this  is  that  a pair  of  hard-way  coupled  concentric 
modules  will  have  a zero  Petzval  sum.  Because  a camera  lens  must 
contain  an  odd  number  of  modules,  the  entire  Petzval  contribution  must 
reside  in  the  final  module.  The  matrix  of  coefficients  of  the  zero- 
Petzval  equation  becomes  singular. 

In  the  modular  design  of  a copy  lens  the  first  and  last  modules 
are  "odd"  modules  in  this  sense  because  all  the  others  are  hard-way 
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couples  and  therefore  contribute  nothing  to  the  Petzval  sum.  It  would 
appear  to  be  possible  to  have  the  Petzval  contribution  of  the  two  ex- 
treme modules  cancel  out  exactly  resulting  in  a concentric  copy  lens 
with  no  field  curvature.  It  turns  out,  however,  that  this  can  work  only 
l i unit  magnification  and  even  then  the  image  distance  becomes  virtual. 
This  feature  certainly  disqualifies  the  design  as  a copy  lens;  it 
might  be  useful  as  a loupe  --  although  even  in  that  case  there  is  no 
particular  advantage  in  having  a zero  field  curvature. 

In  the  concentric  form  of  the  module  all  pupils  coincide  with 
the  common  center  of  curvature.  One  can  ask  if  it  is  not  possible  for 
a similar  event  to  occur  for  some  particular  parameter  value.  A 
module  consists  of  two  spherical  refracting  surfaces.  Suppose  that  an 
aplanatic  point  of  the  first  surface  coincides  with  an  aplanatic  point 
of  the  second.  Then  the  remaining  two  aplanatic  points  would  be  con- 
jugates of  one  another.  If  the  two  module  pupils  coincided  with  these, 
then  we  would  have  a module  configuration  distinct  from  the  concentric 
case  but  with  many  of  its  advantages  and  with  the  possibility  that 
one  could  build  up  a lens  system  for  which  Petzval  could  be  satisfied. 

It  turns  out,  unfortunately,  that  this  configuration  is  inconsistent 
with  the  basic  property  of  a module  --  that  it  be  free  of  spherical 
aberration. 

This  does  suggest  that  it  may  be  possible  to  define  an  entirely 
different  kind  of  module,  one  for  which  third-order  spherical  aberration 
is  not  corrected  but  that  contains  the  aplanatic  configuration  as  a spe- 
cial shape  parameter  value.  This  might  be  important  in  broadening  the 
number  of  building  blocks  available  for  use  in  modular  optical  design. 

In  addition  to  adding  the  "concentric  point"  to  the  list  of 
critical  values  to  be  calculated  in  our  computer  programs  a number  of 
other  more  important  changes  have  been  made.  In  previous  reports  we 
have  described  how,  "in  principle"  distortion  and  lateral  color  could 
be  made  zero  simultaneously  before  the  other  aberrations  were  corrected. 
However,  because  the  two  equations  were  highly  nonlinear  the  task  was 
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by  no  means  an  easy  one.  In  most  cases  we  could  correct  one  or  the 
other  of  these  aberrations,  but  not  both  together.  During  the  past 
quarter  we  have  developed  an  iterative  procedure,  dubbed  RELAX  appro- 
priately enough,  which  does  exactly  this.  Working  on  the  shape  parameters 
of  the  last  two  modules  it  finds  values  for  which  distortion  and  lateral 
color  are  zero. 

This  is  shown  in  Fig.  8,  a print-out  from  program  CAM2  that 
generates  camera  designs  from,  in  this  case,  five  modules.  This  lens 
is  nominally  a 40  in.  f/ 10  camera  lens  with  a semi-field  angle  of  20°. 

The  top  five  lines  (A)  give  the  case,  parameter  value,  and  category  of 
each  of  the  five  modules.  Next  (B)  come  the  coefficients  of  the  power 
parameters  in  the  equations  for  zero  coma,  zero  Petzval,  and  zero  axial 
color.  The  equations  for  distortion  and  lateral  color  do  not  depend 
on  the  power  parameters  and  can  be  calculated  once  and  for  all  at  this 
stage.  The  values  appear  at  (C). 

With  five  modules  any  pair  of  these  three  aberrations  fljay  be  made 
zero  simultaneously.  Here  we  choose  coma  and  Petzval.  The  five  values 
of  the  power  parameters  are  printed  at  (D)  and  the  corresponding  optical 
design,  at  (E),  obtained  from  the  power  parameters  and  the  canonical 
optical  parameters  calculated  elsewhere  in  the  program.  The  program 
stores  these  data  in  the  sequential  file  ZC0PE1  designated  at  (F). 

Because  the  equation  for  zero  coma  is  a quadratic  we  can  expect 
two  distinct  solutions.  A second  is  seen  at  (G) . However,  the  computer 
determines  that  the  design  corresponding  to  this  set  of  power  parameters 
has  a negative  thickness  in  the  first  module,  and  the  design  is  suppressed. 

Note  that  this  program  automatically  finds  the  stop  location. 
Wherever  a pupil  is  real  (rather  than  virtual)  and  in  air,  the  program 
will  print  out  its  location. 

Figure  9 shows  the  printout  of  first  and  third-order  calculations 
from  program  THRAYA.  I believe  it  is  self-explanatory.  Note  that  under 
"sums  of  aberrations"  at  the  bottom  of  the  page  only  axial  c^ior  is 
not  equal  to  zero. 
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Fig.  8,  Printout  from  program  CAM2. 
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Fig.  9.  Printout  of  first-  and  third-order 
calculations  from  program  THRAYA. 
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Time-Sharing  Computer  Program 

At  the  present  time  we  are  working  on  only  one  addition  to  our 
repertoire,  program  11EPTA  for  calculating  third-order  image  and  pupil 
aberration  coefficients,  fifth-order  aberration  coefficients,  and 
seventh-order  spherical  aberration  coefficient.  In  addition,  it  will 
calculate  primary  and  secondary  chromatic  aberration  coefficients. 


Diffraction  Theory 
(A.  S.  Marathay) 


A theoretical  basis  was  formulated  to  study  the  diffraction  of 
light  from  apertures  whose  amplitude  transmittance  is  specified  on  a 
portion  of  a sphere,  The  theory  calculates  the  diffraction  pattern  over 
the  surface  of  an  observation  sphere  concentric  with  the  aperture  sphere. 
This  is  analogus  to  conventional  diffraction  theory  where  the  aperture 
is  situated  in  a plane  and  the  diffraction  pattern  is  calculated  over 
a parallel  plane  at  some  distance.  In  this  formulation  an  appropriate 
Green's  function  is  constructed  by  use  of  the  eigenfunctions  of  the 
Helmholtz  equation  for  spherical  geometry  such  that  it  vanishes  on  the 
sphere.  The  diffracted,  amplitude  is  expressed  as  an  integral  on  the 
surface  of  the  sphere  that  involves  only  the  amplitude  distribution  in 
the  aperture.  In  terms  of  the  Laplace  coefficients  the  integral  re- 
duces to  a product  relation  that  allows  us  to  define  the  transfer 
function  of  free  space.  The  normal  derivative  of  the  Green's  function 
plays  the  role  of  the  impulse  response  of  the  diffraction  problem.  We 
are  going  to  devise  appropriate  computer  programs  to  allow  us  to  cal- 
culate typical  diffraction  patterns  in  this  theory  and  compare  them 
with  the  conventional  diffraction  patterns.1 

In  the  Fourier  space  formulation  of  diffraction  theory  for  aper- 
tures in  a plane  it  is  necessary  to  use  Weyl's  formula.  This  formula 
gives  the  Fourier  transform  of  the  space  part  of  an  expanding  spherical 
wave.  The  derivation  of  Weyl's  formula  is  much  involved.  We  were 
able  to  show  that  it  may  be  derived  as  a natural  consequence  of  the 
inhomogenous  Helmholtz  equation  satisfied  by  the  Green's  function  of 
the  problem.  This  simpler  derivation  gives  more  insight  into  the 
diffraction  problem  and  is  also  pedagogical ly  advantageous. 
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Radiometry 

The  fundamental  laws  of  (noncoherent)  radiometry  have  never 
been  derived  from  the  basic  (wave)  theory  of  light.  Several  attempts 
have  been  made  in  the  literature,  notably  by  Walther2  and  Marchand  and 
Wolf,3  to  derive  these  laws.  Their  formulations  are  unsatisfactory 
for  several  reasons: ‘ (1)  they  introduce  certain  mathematical  functions 
without  any  physically  measurable  properties--they  are  merely  mathe- 
matical constructs,  and  (2)  their  results  indicate  that  for  a nonco- 
herent source  the  radiant  intensity  behaves  as  cos20  where  0 is  the 
angle  that  the  detector  makes  with  the  normal  to  the  surface.  That  is, 
they  show  that  a noncoherent  source  is  not  Lambertian.  This  state  of 
affairs  is  unacceptable. 

We  therefore  felt  it  was  necessary  to  reformulate  the  theory  to 
generalize  the  laws  of  radiometry  to  partially  coherent  fields,  such 
that  it  will  reproduce  the  well-known  laws  of  radiometry.  We  present 
a new  approach  to  generalizing  the  laws  of  radiometry  to  partially 
coherent  fields.  A new  complex  cross-correlation  function  F is  in- 
troduced that  depends  on  two  space  points,  two  directions,  and  the 
time  difference  t.  This  function,  F,  and  the  mutual  coherence  function, 
r,  are  the  basis  for  defining  the  generalized  radiometric  quantities. 
Special  forms  of  the  generalized  radiance  L are  studied.  These  defini- 
tions are  applied  to  two  special  cases.  In  particular  we  show  that  a 
uniform,  noncoherent  source  of  infinite  extent  is  Lambertian.  In  partic- 
ular for  this  source  our  generalized  functions  reduce  to  the  traditional 
radiometric  quantities  and  we  also  successfully  rederive  their  inter- 
relationships. The  relationships  are4 

♦o  = 

M o E * L o 
and 

Jo  = jf  L$  cos0, 


where  $ is  the  flux  with  units  of  W,  Mq  is  the  radiant  emittance 
W m"2,  Jo  is  the  radiant  intensity  W sr"1,  and  the  radiance  L o has  the 
units  of  W m"2  sr"1.  The  quantity  ji  is  the  area  element  of  the  source. 
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Image  Plane  Speckle 
(M.  J.  Lahart  and  A.  S.  Marathay) 

We  have  finished  the  main  theoretical  aspect  of  our  work  con- 
cerning reduction  of  speckle  by  using  weak  diffusers.  Our  calculations 
have  been  based  on  the  study  of  second-order  statistics.  In  particular, 
we  calculated  the  intensity  correlation  function  of  the  speckle  pat- 
tern as  a function  of  the  aperture  size  of  the  imaging  system  and  of  the 
characteristics  of  the  diffuser. 

Diffusers  are  used  in  coherent  optical  systems  to  reduce  dif- 
fraction noise  from  dust  or  blemishes.  They  cause  the  dust  particle 
or  noise  source  to  be  illuminated  from  several  directions  so  that  the 
effects  of  noise  are  spread  out  or  averaged  and  generally  made  less 
apparent.  The  process  of  averaging  can  create  speckle  or  some  other 
artifact,  however,  and  designers  of  coherent  systems  must  know  how  to 
control  both  speckle  and  diffraction  noise  and  reduce  their  combined 
effect. 

We  have  investigated  the  way  diffraction  noise  is  reduced  by  a 
diffuser  by  computing  its  effect  on  several  individual  noise  patterns 
produced  by  small  particles.  Noise  of  this  kind  is  nearly  eliminated  if 
light  is  scattered  through  a relatively  small  minimum  angle,  and  it  is 
reduced  only  li<  le  nu-'-e  if  light  is  scattered  through  much  larger 
angles.  The  coritr.  v of  the  speckle  pattern  may  be  high  if  light  is 
scattered  through  la  ;er  angles  than  necessary,  however,  unless  the 
aperture  size  is  also  large.  The  minimum  scattering  angle  is  an  impor- 
tant design  parameter  because  it  determines  how  large  the  system  aper- 
ture must  be  to  restrict  speckle  contrast  to  an  acceptable  level. 

Recently  it  has  been  i.  cognized  that  random  phase  plates  can  be 
used  to  scatter  light  through  small  angles.  The  speckle  pattern  that 
is  produced  is  the  spatially  filtered  image  of  the  phase  object,  and 
its  contrast  is  small  if  nearly  all  of  the  light  is  transmitted  by  the 


system.  Phase  plate  diffusers  can  permit  light  to  be  scattered  through 
the  minimum  angle  necessary  to  reduce  diffraction  noise,  without  scattering 
it  through  the  larger  angles  that  make  speckle  more  apparent. 

A random  phase  plate  scatters  light  through  large  angles  like  a 
rtrong,  ground  glass  diffuser  if  the  phase  variance  is  large  and  the 
phase  correlation  distance  is  small.  Light  may  be  scattered  through 
smaller  angles  if  the  phase  correlation  length  is  large,  but  the  scattering 
distribution  may  be  coarse  and  uneven,  especially  if  the  correlation 
length  is  comparable  to  the  dimensions  of  the  diffuser.  If  the  phase 
correlation  distance  is  small,  the  phase  variance  must  also  be  small. 
However,  if  the  variance  is  too  small,  the  specular  component  of  the 
radiation,  the  light  that  is  transmitted  without  being  scattered,  may  be 
large,  so  that  the  effectiveness  of  the  diffuser  is  reduced.  A phase 
plate  that  is  used  as  a noise  reducing  diffuser  should  have  r.  small 
specular  component,  but  a small  enough  phase  variance  so  tha’  the 
correlation  distance  can  be  small. 

Oi.i  work  has  shown  that  if  the  phase  of  a random  diffuser  is 
normally  distributed,  the  percentage  of  specularly  transmitted  light  is 
small  when  the  phase  standard  deviation  is  one  or  two  adians.  Not  only 
is  almost  all  of  the  light  scattered  in  this  case,  but  the  detailed 
characteristics  of  th^  scattering  distribution  can  be  specified  because 
the  distribution  is  nearly  proportional  to  the  phase  Wiener  spectrum. 

These  small  phase  standard  deviations  appear  to  be  the  optimum  for 
noise  reducing  diffusers  in  imaging  systems. 

We  have  investigated  the  image  speckle  pattern  of  weak  diffusers 
in  detail.  The  speckle  contrast  can  be  surprisingly  high  even  when 
most  of  the  diffracted  light  is  transmitted,  but  it  sometimes  consists 
primarily  of  high  spatial  frequencies.  (This  is  not  true  of  speckle 
patterns  of  strong  diffusers.)  The  optical  system  can  sometimes  be 
designed  to  eliminate  speckle  noise  that  is  not  in  the  same  part  of  the 
spatial  frequency  spectrum  as  the  image;  the  speckle  may  not  be  objec- 
tionable when  this  is  the  case. 
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The  statistical  law  that  the  phase  obeys  was  found  to  have  an 
important  influence  on  the  speckle  pattern.  Diffusers  whose  phase  was 
negative  exponentially  distributed  created  higher  contrast  patterns 
and  scattered  less  light  than  those  whose  phase  was  normally  distributed. 
We  found  that  techniques  for  making  negative  exponential  diffusers  can 
be  modified  to  produce  normally  distributed  diffusers  with  prescribed 
phase  Wiener  spectra. 

We  investigated  the  influence  of  the  coherence  of  the  illumination 
on  the  speckle  pattern.  The  speckle  contrast  is  small  if  the  light  co- 
herence length  is  sufficiently  short,  but  moderately  incoherent  light 
can  produce  a speckle  pattern  whose  Wiener  spectrum  is  dominated  by  low 
spatial  frequencies  but  whose  contrast  is  nearly  the  same  as  a coherent 
lignt  pattern.  Partially  coherent  illumination  is  often  proposed  as  a 
solution  to  the  speckle  problem;  our  work  shows  that  its  use  is  not 
always  a good  solution,  and  it  can  even  make  the  problem  worse. 

Some  work  was  necessary  to  delineate  the  regions  of  validity  of 
our  computations  because  most  of  our  calculations  are  ensemble  averages. 
The  ergodic  theorem  guarantees  that  a single  pattern  can  be  characterized 
by  its  ensemble  average  if  it  is  large  enough,  but  this  characteri- 
zation may  not  be  appropriate  if  the  pattern  is  of  finite  size.  We 
analyzed  the  statistical  accuracy  of  finite  pattern  calculations  and 
showed  that  ensemble  averages  are  accurate  under  many  but  not  all 
practical  conditions,  and  we  established  minimum  sizes  of  formats  under 
which  they  are  good  descriptors. 


OPTICAL  TRACKING 
(J.  C.  Wyant) 


Optical  Tracking  System 

The  optical  tracking  system  study  for  the  Space  and  Missile 
Test  Center  (SAMTEC)  at  Vandenberg  Air  Force  Base  has  continued  to  ful- 
fill the  seven  technical  tasks  described  in  the  past  SAMSO  quarterly 
report.  An  interim  report  was  written  giving  the  progress  to  date.  The 
report  describes  recent  work  in  optical  detection  systems  that  establishes 
the  limits  of  accuracy  and  performance  capability  of  state-of-the-art 
optical  tracking  systems  for  collecting  metric  data  on  the  Western  Test 
Range.  While  imaging  systems  cannot  provide  the  desired  metric  data, 
they  can  aid  in  the  tracking  of  a laser  radar  system  and  can  also  give 
engineering  data.  As  long  as  the  missile  is  equipped  with  a retro- 
reflector,  a laser  ranger  can  give  the  missile  range  to  within  a 10-cm 
error,  and  a laser  Doppler  system  can  give  the  missile  velocity  to  within 
an  error  of  0.01  ft/sec  or  less  during  the  bum  stages  of  the  missile 
any  time  of  the  day  or  night,  as  long  as  the  sun  is  not  directly  within 
the  field  of  view  and  the  weather  is  clear. 

Use  of  White-Light  Extended  Source  with  a Lateral  Shear 
Interferometer  for  Real-Time  Wavefront  Correction  Systems 

An  analysis  was  performed  to  determine  the  accuracy  with  which  an 
ac  heterodyne  lateral  shear  interferometer  can  measure  wavefront  aber- 
rations if  a white-light  extended  source  is  used  with  the  interferometer 
and  shot  noise  is  the  predominant  noise  source.  The  analysis  shows  that 
for  uniform  circular  or  square  sources  larger  than  a derived  minimum 
size,  the  wavefront  measurement  accuracy  depends  on  only  the  radiance 
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of  the  source  and  not  on  the  angular  subtense  of  the  source.  For  a 
radiance  of  10  W/m2  sr  the  rms  wavefront  error  is  approximately  1/30  wave, 
assuming  the  singal  is  shot  noise  limited.  For  both  uniform  circular 
and  square  sources  an  optimum  shear  distance  is  approximately  half  the 
aperture  diameter  required  to  resolve  the  light  source.  The  analysis  is 
of  particular  importance  because  of  the  use  of  lateral  shear  interferom- 
eters in  active  optics  systems,  especially  for  systems  in  which  wavefront 
deformations  introduced  by  atmospheric  turbulence  are  being  corrected. 

Etched  Gratings 

Work  is  continuing  on  the  fabrication  of  etched  gratings  in  glass. 
The  techniques  are  believed  useful  for  active  optics  systems  in  which 
changes  in  the  figure  of  the  primary  mirror  are  being  measured.  Because 
the  gratings  are  etched  in  the  mirror  surface  directly,  there  is  no 
problem  in  the  grating  moving  relative  to  the  mirror  surface.  Because 
the  etch  depth  is  small  and  the  diffraction  orders  are  out  of  the  field 
of  view,  the  etched  grating  will  not  affect  the  performance  of  the  optical 
system.  Progress  during  the  past  quarter  includes  improving  the  repeat- 
ability of  the  grating-forming  process  and  improving  the  uniformity 
of  a 2-in. -square  grating.  Work  is  presently  being  performed  to  make  a 
6-in.-diam  grating  on  a parabolic  mirror. 

Comparison  of  Heterodyne  and  Direct  Detection  Systems 

An  analysis  is  being  performed  to  compare  the  signal-to-noise 
ratios  of  heterodyne  and  direct  detection  systems  in  the  presence  of 
atmospheric  turbulence  and  several  background  radiation  levels.  The 
particular  wavelength  of  interest  is  10.6  urn,  although  the  use  of  shorter 
wavelengths  is  also  being  investigated. 
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QUANTA!  OPTICS 

Picosecond  Spectroscopy 
(J.  C.  Matter,  M.  0.  Scully,  A.  L.  Smirl) 

The  measurement  of  the  intraband  lifetime  of  carriers  in  semi- 
conductors has  shifted  from  germanium1  to  mercury  cadmium  telluride, 

Hg.  Cd  Te.  This  material  was  chosen  chiefly  for  two  reasons: 

(1)  HgCdTe  is  a direct  band  gap  semiconductor2  in  contrast  to  indirect 
gap  germanium,  which  may  provide  a simpler  model  for  theoretical  analysis; 
and  (2)  the  width  of  the  band  gap  can  be  selected  by  varying  the  alloy 
parameter  x,3  the  mole  fraction  of  cadmium  telluride  (CdTe) . 

We  do  not  have  definitive  experimental  results  to  present  at 
this  time.  It  has  become  apparent  that  surface  preparation,  i.e., 
etching,  is  an  important  factor  for  HgCdTe,  while  for  germanium  it  is 
not  (for  these  experiments).  This  seems  reasonable  in  light  of  the 
fact  that  HgCdTe  is  a much  softer  solid  than  Ge.  Thus,  when  the  sample 
is  being  polished,  surface  damage  may  be  much  more  severe  for  HgCdTe 
than  for  Ge. 

We  will  briefly  describe  the  experiments  performed  thus  far. 

For  unetched  samples  of  HgCdTe,  no  saturation  of  transmission  at  high 
intensities  (X  = 1.06  pm)  was  observed  and  subsequently  no  lifetime  was 
measured.  For  HgCdTe  samples  etched  on  the  front  (incident  radiation) 
surface,  saturation  of  transmission  was  observed  for  one  sample  of  large 
x value.  Lifetime  curves  were  obtained  as  a function  of  temperature 
and  as  a function  of  excitation  pulse  intensity.  In  order  to  determine 
if  this  observed  lifetime  is  dominated  by  the  unetched  back  (mounted 
on  KZF-2  glass  substrate)  surface,  we  are  now  preparing  new  samples 
that  will  be  etched  on  both  surfaces. 
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Transient  Interactions 


(R.  L.  Shoemaker  and  E.  Van  Stryland) 

During  the  past  quarter  we  have  successfully  made  direct 
measurements  of  transition  dipole  matrix  elements  in  a molecule  using 
optical  nutation.  The  transition  dipole  matrix  element,  vah  - \ \ erpyh , 
is  an  important  parameter  in  any  spectroscopic  experiment  since  it  is 
In^l2  that  fundamentally  determines  the  strength  of  a transition  from 
the  state  ip  to  the  state  \p^.  Accurate  values  of  are  useful  in  a 
wide  variety  of  applications  including  the  study  of  population  dis- 
tributions in  gases,  radiative  energy  transfer  problems,  and  the  search 
for  new  laser  transitions.  In  addition,  for  vibrational  transitions  in 
molecules,  the  vibrational  piece  of  can  be  used  as  a test  for  vi- 
brational potentials  and  wave  functions.  Thus  accurate  methods  for 
measuring  transition  dipole  matrix  elements  are  of  considerable  interest. 

We  have  demonstrated  how  direct,  accurate  measurements  of  tran- 
sition dipole  matrix  elements  can  be  made  using  optical  nutation.  This 
phenomenon  is  seen  as  a damped  oscillation  in  the  transmitted  intensity 
that  occurs  when  a molecular  transition  is  suddenly  brought  into  re- 
sonance with  a laser  beam.  The  oscillation  arises  because  molecules  are 
driven  coherently  back  and  forth  between  the  ground  and  excited  states 
producing  an  alternating  absorption  and  emission  of  radiation.  If  the 
transition  is  excited  exactly  on  resonance  by  a uniform  optical  ield 
with  amplitude  Eq,  the  oscillation  frequency  is  just  Thus  a 

measurement  of  the  oscillation  frequency  together  with  a measurement  of 
the  laser  intensity  I * '-joeo^o2  gives  the  matrix  element  directly. 

Transition  dipole  matrix  elements  can  also  be  determined  from 
curve  of  growth  measurements  where  the  absolute  absorption  as  a function 
of  pressure  is  measured.  This  method  can  give  accurate  results,  but 
only  when  certain  conditions  are  met.  Optical  nutation  is  complementary 
to  the  curve  of  growth  method  in  the  sense  that  the  accuracy  of  the  two 


techniques  depends  on  different  parameters  so  that  one  can  often  be  used 
in  situations  where  the  other  fails. 

Our  optical  nutation  measurements  were  made  on  the  deuterated 
ammonia,  NH2D,  using  a Stark  switching  technique  described  below.  Tran- 
sition dipole  matrix  elements  are  obtained  for  two  vibration-rotation 
transitions  in  the  molecule,  and  since  the  transition  assignments  are 
known,  the  vibrational  piece  of  the  transition  dipole  matrix  elements 
along  two  different  axes  in  the  molecule  can  be  determined.  These  quan- 
tities should  provide  a sensitive  tost  of  the  V2  vibrational  wave  function, 
in  NH2D. 

There  has  been  one  previous  attempt  to  measure  a transition  dipole 
matrix  element  using  optical  nutation.  Hocker  and  Tang  made  a mea- 
surement on  SF6  using  a pulsed  laser.  It  is  hard  to  assess  the  accuracy 
of  their  results  because  of  the  difficulties  inherent  in  a pulsed  laser 
experiment.  In  particular,  the  nutation  frequency  varies  because  the 
optical  field  strength  changes  during  the  pulse  and  the  laser  frequency 
can  be  chirped.  It  is  also  difficult  to  measure  accurately  the  nutation 
signal  that  does  occur,  since  it  appears  as  a small  signal  on  a time 
varying  background.  Furthermore,  a determination  of  the  optical  field 
strength  is  extremely  difficult  since  the  pulse  shape,  the  pulse  energy, 
and  the  transverse  profile  of  the  beam  must  all  be  reproducible  and 
accurately  known.  Thus  the  measurement  of  transition  dipole  matrix 
elements  using  this  approach  does  not  appear  very  promising  if  accuracies 
much  beyond  an  order  of  magnitude  estimate  are  desired. 

The  difficulties  just  mentioned  can  be  avoided,  however,  by  using 
the  Stark  switching  technique  developed  by  Brewer  and  Shoemaker.  In  this 
method,  which  we  have  used  for  the  measurements  reported  here,  a cw  laser 
is  utilized  and  molecular  transitions  are  shined  in  and  out  of  resonance 
using  the  Stark  effect.  Thus,  one  merely  passes  the  laser  beam  through 
the  sample  and  applies  a pulsed  electric  field  across  it.  The  optical 
nutation  signal  is  observed  by  m.ui.’toring  the  transmitted  light  with  a 
photodetector.  Using  a pulsed  electric  field  instead  of  a pulsed  laser 
produces  several  significant  advantages.  First,  because  the  laser  beam 


is  cw,  it  produces  only  a dc  signal  in  ♦'he  detector.  The  only  ac  signal 
present  is  just  the  optical  nutation  signal  one  wants  to  observe!  This 
allows  the  nutation  signal  to  be  accurately  measured.  Second,  since 
the  pulsed  electric  field  can  be  a nearly  ideal  step  function  and  the 
cw  laser  is  stable  in  amplitude  and  frequency,  an  exact  theoretical  cal- 
culation of  the  optical  nutation  signal  can  be  made  with  p^  as  the  only 
adjustable  parameter.  Finally,  the  laser  intensity  may  be  accurately 
measured  by  using  a calibrated  thermopile  to  obtain  the  total  beam  power 
and  a scanning  aperture  to  obtain  the  beam's  transverse  intensity  profile. 
Thus  one  expects  that  accurate  measurements  of  p^  can  be  made,  and  we 
have  found  that  these  expectations  are  fulfilled. 

A fundamental  limitation  of  the  Stark  switching  technique  is,  of 
course,  that  transitions  with  a substantial  Stark  effect  must  be  used. 
However,  the  technique  can  be  generalized  by  frequency  switching  the  laser 
instead  of  the  sample.  Optical  nutation  can  then  be  observed  in  any 
transition  that  overlaps  a laser  line. 

We  have  developed  a semiclassical  treatment  of  optical  nutation 
that  can  be  quantitatively  compared  with  our  experiments.  We  wish  to 
consider  what  happens  when  a cw  laser  beam  passes  through  a low  pressure 
gas  sample  with  a resonant,  nondegenerate  transition.  The  laser  fre- 
quency is  assumed  initially  to  be  somewhere  within  the  Doppler 
broadened  linewidth  of  the  transition  (-  80  MHz  FWW1  for  NH2D  at  X ■ 10  pm). 
Only  those  molecules  with  a z component  of  velocity  v ',  such  that  they 
are  Doppler  shifted  into  resonance,  are  excited,  and  this  narrow  velocity 
group  has  a homogeneous  linewidth  that  is  typically  100  kHz  tc  1 MHz  FWHM. 

At  t * 0,  a step  function  electric  field  is  applied  acrosr  the  sample 
that  shifts  the  center  frequency  of  the  transition  to  a new  value.  If  the 
frequency  shift  is  several  homogeneous  linewidths  or  more,  a new  velocity 

group,  v ",  will  suddenly  be  shifted  into  resonance  producing  an  optical 

z 

nutation  signal.  In  addition  to  the  optical  nutation,  there  will  also  be 
an  optical  free  induction  decay  signal  (FID)  from  the  velocity  group,  i>2', 
that  was  shifted  out  of  resonance.  This  signal  can  be  ignored  in  practice 
because  the  FID  appears  as  a heterodyne  beat  with  the  laser  and  decays 
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as  cxpt-Y*  - |/Y2  * {vajE^/’h)2t'\,  where  y is  the  collisional  decay  rate. 
Thus  the  FID  has  essentially  died  out  by  the  end  of  the  first  cycle  c f 
optical  nutation,  and  even  during  that  cycle,  its  amplitude  can  be  made 
negligibly  small  by  simply  Stark  shifting  far  enough  to  make  the  heter- 
odyne beat  frequency  lie  outside  the  bandwidth  of  the  detection  system. 

In  light  of  this  we  may  neglect  the  previously  excited  velocity  group, 

' , and  consider  only  the  response  of  an  unexcited  transition  suddenly 
shifted  into  resonance  at  t = 0.  Our  approach  is  to  calculate  the  molec- 
ular polarization  by  means  of  the  density  matrix  equations  of  motion, 
and  then  to  use  Maxwell's  equation  to  find  the  signal  that  appears  at 
the  detector. 

We  take  the  optical  field  to  be 

E^  * Eq  cos(n^t  - kz) . (1) 

Considering  now  an  ensemble  of  molecules  with  axial  velocities  between 
v2  and  v * do  , and  making  a (nonrelativistic)  transformation  to  a 
moving  coordinate  frame  in  which  these  molecules  are  at  rest,  this 
becomes 


* E0  cosC^'t  - kz')  f 


(2) 


where  - kv g and  z'  is  the  molecular  position  in  the  moving 

frame.  We  assume  an  optically  thin  sample  so  that  Eq  does  not  depend 
on  z.  This  condition  is  certainly  satisfied  in  our  experiments  chore  the 
peak  absorption  is  -0.3%. 

The  equation  of  motion  for  the  molecular  density  matrix  can  be 

written 


ir  1,-+  -*-0. 

-— (p  - P ) 


(3) 


where  p is  the  density  matrix  in  thermal  equilibrium  and  we  have  intro- 
duced a single  phenomenological  decay  time  t.  This  is  equivalent  to  the 
assumption  that  the  molecular  relaxation  processes  are  dominated  by 
strong  collisions  in  which  the  collisional  transition  probability  to  any 


level  is  simply  proportional  to  the  Boltzmann  factor  for  that  level. 
Such  an  assumption  is  not  in  fact  justified  for  the  NH2D  transitions 
studied  here  since  the  relaxation  time  for  the  population  difference, 
f>bb"paa‘  an<*  t*ie  relaxat*on  t*me  ^or  off-diagonal  element 
differ.  However,  for  low  sample  pressures  the  decay  of  the  observed 
optical  nutation  signal  is  dominated  by  effects  due  to  the  transverse 
intensity  profile  of  the  laser  beam.  Thus  the  collisional  decay  model 
one  uses  does  not  significantly  affect  the  results  of  the  calculation, 
and  we  can  proceed  using  the  simple  strong  collision  model. 

Upon  solving  Eq.  (3),  calculating  the  resulting  polarization, 
and  using  Maxwell's  equations  to  find  the  transmitted  intensity,  we 
obtain 
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where  Jq  * cz.qEqz/2  and  g = ^[(fcv^)2  ♦ In  the  limit  as 

t •*  ®,  this  equation  reduces  to  the  standard  expression  for  steady- 
state  saturated  absorption  at  the  peak  of  a Doppler  broadened  line.  For 
short  times  exp(-t/-r)  -*•  1,  and  the  integral  in  Eq.  (4)  can  be  evaluated 
approximately  to  give  2t<7q[  (u^Eo/bjt].  Since  the  zero  order  Bessel 
function  Jq  is  approximately  unity  for  small  t,  the  initial  peak  value 
of  the  optical  nutation  signal  will  be 
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We  note  that  Eq.  (5)  is  exactly  the  expression  one  obtains  for  uneaturated. 
absorption  assuming  a small  optical  field.  This  is  not  unexpected  since 
for  a short  time  after  the  molecules  are  switched  into  resonance,  the 
level  populations  have  not  yet  changed  appreciably  and  one  sees  the  full 
unsaturated  absorption. 
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Experiment 

A block  diagram  of  the  experimental  arrangement  is  shown  in 
Fig.  10.  A 1.5  m grating-controlled  CO2  laser  is  used  that  produces 
approximately  2 to  4 W,  TEMqo  mode,  on  any  one  of  approximately  50  C02 
lines.  A small  650  Hz  dither  is  applied  to  the  output  mirror  producing 
a variation  in  the  discharge  tube  impedance.  This  signal  is  used  to 
lock  the  laser  to  the  peak  of  its  gain  curve  with  a long-term  frequency 
stability  of  ±1  MHz.  The  short-term  frequency  jitter  is  less  than  100  KHz. 
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Fig.  10.  Block  diagram  of  the  experimental  arrangement. 
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The  NH2D  gas  samples  are  held  in  a Stark  cell  equipped  with 
NaCl  windows  and  containing  two  4 x 30  cm  long  glass  plates  metallized 
on  their  inner  surfaces  and  separated  by  .5946  t .0005  ram  quartz  spacers. 
After  passing  between  the  Stark  plates,  the  transmitted  laser  beam  is 
monitored  by  Au:Ge  photoconductive  detector,  that  in  combination  with  a 
preamplifier,  has  a bandwidth  of  approximately  10  MHz.  A pulse  generator 
supplies  step  function  voltage  pulses  across  the  Stark  plates  at  a 
repetition  rate  of  approximately  20  kHz,  and  the  detector-preamp  output  is 
fed  into  a PAR  162  Boxcar  integrator  that  provides  signal  averaging 
and  displays  the  results  on  an  x-y  recorder. 
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NH2D  was  synthesized  by  simply  mixing  NH3  and  D2O.  The  H and 
D atoms  in  both  NH3  and  D2O  exchange  rapidly  among  themselves,  resulting 
in  a statistical  mixture  of  partially  deuterated  water  and  ammonia. 

This  mixture  was  vacuum  distilled  several  times  holding  the  mixture  at 
-78°C  and  collecting  the  distillate  at  -196°C.  Since  the  vapor  pressure 
of  H2O  is  10" 3 Torr  at  -78°C,  as  compared  to  50  Torr  for  NH3,  all 
water  can  be  removed  in  this  manner.  The  purified  ammonia  mixture 
contains  45%  NH2D  if  the  NH3:D20  ratio  is  1.33:1  by  volume. 

We  first  examined  the  ) - (1,  4o  4,  4)  + (2,  5x4,  5)  NH2D 
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transition  that  was  Stark  tuned  approximately  2300  MHz  into  resonance 
with  the  P(14)  CO2  laser  line  at  949.48  cm-1.  The  Stark  field  is  per- 
pendicular to  the  optical  field  direction  such  that  AMj  = ±1  selection 
rules  apply,  the  bMj  * ±1  transition  being  by  far  the  more  intense. 

This  could  readily  be  made  more  accurate  in  future  experiments.  We 
also  measured  y^  using  a lower  laser  power,  1.28  W at  3.9  mTorr. 

Here  the  signal  is  correspondingly  more  sensitive  to  x because  the 
nutation  frequency  is  lower.  However,  the  fit  is  still  good  (see 
Fig.  11)  and  yielded  a value  of  y^  * .0356  t .0021  Debye.  This  differs 
from  the  higher  power  measurement  by  only  2%.  The  average  value  for 
the  dipole  moment  is  y^  * .0353  t .0021  Debye.  We  also  made  measurements 
on  two  M components  of  the  (v,<7)  = (0,4o4>  *►  (l,5os)  transition  that 
may  be  Stark  tuned  into  resonance  with  the  P(20)  CO2  laser  line  at 
949.19  cm'1.  This  transition  involves  a three- level  system  arising  from 
a small  Stark  and  quadrupole  splitting  of  the  upper  level.  This  splitting 
causes  a quantum  beat  in  the  nutation  signal.  As  discussed  in  a previous 
report,  we  have  calculated  the  expected  nutation  signal  in  this  three- 
level  system  and  found  the  possible  deviatio)  in  the  value  of  y^  from 
the  two-level  theory  to  be  less  than  15%.  Tho  deviation  should  be 
small  since  one  of  the  upper  states  has  a much  smaller  dipole  matrix 
element  connecting  it  to  the  lower  level.  The  values  obtained  for  y^ 
are 

y^  ■ .0355  for  Mj  ■»  4 + 5,  and 

uab  * 


.0295  for  Mj 


3+4. 
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Taking  into  account  the  Mj  value  degeneracy  factors  these  values  agree 
to  within  7%. 


TIME  [p sec) 

Fig.  11.  Experimental  nutation  signal  using  3.9  mTorr  of 
NH2D  and  a laser  output  power  of  1.28  W. 


From  these  /alues  for  the  total  dipole  moments,  the  vibrational 
piece  is  extracted  as  follows:  the  upper  and  lower  state  wave  functions 

are  expanded  in  terms  of  the  symmetric  rotor  basis  wave  functions,  and 
uv  is  obtained  from 


W 


v 


^total 

<< b,  > 

lower ' 1 r upper 


where  the  denominator  is  just  the  rotational  direction  cosine  matrix 
element. 


This  transition  involves  a single  nondegenerate  two- level  system 
so  that  the  above-described  theory  applies.  The  experimental  nutation 
signal  using  6.6  mTorr  of  NH2D  and  a laser  output  power  of  2.34  W is 
shown  in  Fig.  12  along  with  the  theoretical  fit.  Note  that  the  detector 
sees  the  field  emitted  by  all  the  molecules  interacting  with  the  laser 
beam  because  of  diffraction  at  10.6  urn.  We  verified  this  by  varying 
the  aperture  directly  behind  the  cell  and  noted  that  the  nutation  signal 
remained  unaltered.  Therefore,  since  the  laser  field  £0.  is  a function 
of  R,  the  r'dial  distance  from  the  center  of  the  beam,  the  signal  S must 
be  the  integral  of  Eq.  (4)  over  R 


S « 2ir  [ I(t)R  dR.  (6) 

J0 

where  I{t)  is  now  a function  of  R.  The  beam  profile,  Eq2(R),  was  very 
nearly  Gaussian  and  was  obtained  by  measuring  the  laser  power  both  in 
front  of  and  behind  the  Stark  cell  using  a small  scanning  aperture  in 
the  vertical  and  horizontal  directions.  These  profiles  were  averaged, 
and  Eq.  (6)  was  numerically  integrated  to  obtain  the  curve  of  Fig. 12. 
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Fig.  12.  Experimental  nutation  signal  using  6.6  mTorr  of 
NH  D and  a laser  output  power  of  2.34  W. 
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The  value  of  t,  the  lifetime,  and  v^EqUi,  where  Eq  is  now  the 
peak  field  intensity,  are  t * 1.2  usee  and  v^Eq/H  ■ 9.0  ± .5  MHz.  The 
errors  are  set  by  a 2%  fitting  error  and  a 5%  beam  area  measure  ,nt  error. 
The  fit  was  very  insensitive  to  the  value  chosen  for  t.  This  is  due  to 
the  fact  that  the  nutation  signal  from  the  molecules  in  the  center  of  the 
beam  (high  power  and  thus  high  frequency)  and  the  edge  of  the  beam  (low 
power  and  low  frequency)  interfere  to  give  a short  lifetime  for  the  nuta- 
tion signal.  From  the  value  of  v^Eq/H  the  dipole  moment  can  be  deter- 
mined if  we  know  Eq.  Eq  is  given  by  the  absolute  power  measurement. 

This  power  was  measured  with  a Coherent  Radiation  model  210  power  meter 
that  had  been  calibrated  to  within  5t  with  a C02  laser  against  a standard 
traceable  to  the  National  Bureau  of  Standards.  We  used  the  average  of 
two  such  meters  that  checked  to  within  4%. 

The  value  of  the  dipole  moment  that  resulted  is  v ^ * .0349  ± .002 
Debye.  The  bulk  of  the  error  here  is  due  to  the  beam  profile  measurement. 

The  value  of  uv  obtained  for  the  P(14)  transition  was  .156  ± .014 
Debye,  whereas  for  the  P(20)  transition,  ■■  ,09V  ± .015.  These  results 
are  the  vibrational  transition  dipole  matrix  elements  along  the  a and  a 
principal  axes  in  the  molecule,  respectively.  To  our  knowledge,  these 
results  represent  the  first  measurements  of  vibrational  dipole  matrix 
elements  ever  made  in  an  asymmetric  rotor  molecule. 
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Electromagnetic  Anisotropy 
(J.  Small) 

An  experimental  test  tor  anisotropy  in  the  velocity  of  light  has 
been  designed  after  the  method  of  E.  W.  Silvertooth. 1 The  experiment 
consists  of  propagating  highly  stabilized  electromagnetic  frequencies 
between  two  points.  The  entire  experiment  is  then  rotated  to  detect 
differential  changes  in  propagation  velocity  with  respect  to  direction. 

A detailed  analysis  shows  that  the  absolute  frequencies  used  are  unimpor- 
tant, but  the  fractional  frequency  stability  of  the  oscillators  contri- 
butes directly  to  the  signal- to- noise  ratio  of  the  experiment. 

Frequency  stability  ic  a precisely  defined  mathematical  concept2 
that  depends  upon  the  sample  time  during  a frequency  measurement.  Dif- 
ferences between  short-term  and  long-term  stability  are  thus  rigorously 
quantified.  In  the  present  test,  the  sample  times  of  interest  are  deter- 
mined by  the  rotation  rate  of  the  experimental  apparatus.  It  is  highly 
important  that  the  oscillators  used  have  adequate  frequency  stability 
during  a single  rotation  period. 

The  oscillators  chosen  for  the  experiment  are  methane-stabilized 
helium-neon  lasers  operating  at  a wavelength  of  3.39  vim.  They  have 
excellent  stability.  "The  frequency  stability  from  1 msec  out  to  ten 
seconds  matches  or  exceeds  the  stability  of  other  oscillators  which 
have  been  used  as  standard  frequency  sources,  e.g.,  hydrogen  masers, 
cesium  beam  devices,  crystal  oscillators,  etc."3  In  addition,  the  short 
infrared  wavelength  eliminates  many  experimental  difficulties,  such  as 
multipath  propagation  and  diffraction  effects,  that  would  be  encountered 
with  microwave  or  radio  frequency  oscillators. 

The  flow  diagram  of  Fig.  13  shows  the  construction  plan  for  the 
experimental  apparatus.  Progress  to  date  is  indicated  by  heavy  under- 
lining. The  work  divides  naturally  into  three  general  areas:  (1)  con- 
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struction  of  the  stable  lasers,  (2)  construction  of  an  air-cushioned 
rotating  table,  and  (3)  computer- control led  data  analysis. 

The  laser  frequency  stabilizing  electronics  are  being  constructed 
§ under  the  supervision  of  the  National  Bureau  of  Standards  in  Boulder, 

Colorado.  When  integrated  with  the  ULE  quartz  laser  cavities  presently 
under  construction  at  the  Optical  Sciences  Center,  f.-equency  stabilities 
equal  to  the  present  state  of  the  art  may  be  expected, 
ft 
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Atomic  and  Molecular  Physics  Laboratory 
(W.  Wing) 

Work  has  continued  on  establishing  the  Atonic  and  Molecular 
Physics  Laboratory  in  the  Physics  Department.  The  equipment  for  three 
experiments  brought  from  Yale  University  in  the  summer  of  1974  is 
being  set  up,  and  reassembly  of  one  experiment  (on  infrared  resonance 
in  molecular  ions)  is  finished.  A dye  laser  pumped  by  a nitrogen 
laser  is  being  constructed  to  excite  Rydberg  states  in  helium  atoms. 

A new  experiment  has  begun  on  the  pumping,  mechanism  of  the  helium- 
neon  laser. 
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